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We consider the evolution and dissipation of vortex rings in a condensate at non-zero temperature, 
in the context of the classical field approximation, based on the defocusing nonlinear Schrodingcr 
equation. The temperature in such a system is fully determined by the total number density and 
the number density of the condensate. A vortex ring is introduced into a condensate in a state 
of thermal equilibrium, and interacts with non-condensed particles. These interactions lead to a 
gradual decrease in the vortex line density, until the vortex ring completely disappears. We show 
that the square of the vortex line length changes linearly with time, and obtain the corresponding 
universal decay law. We relate this to mutual friction coefficients in the fundamental equation of 
vortex motion in superfluids. 

PACS numbers: 03.65.Sq, 03.75.Kk, 05.65.+b, 67.40.Vs, 67.57.De 
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The processes of self-organization, formation of large- 
scale coherent localized structures and interactions of 
these structures with small-scale fluctuations are at the 
heart of nonlinear sciences, ranging from classical tur- 
bulence, superfluids, ultracold gases and Bose-Einstein 
Condensates (BECs) , to the formation of the early Uni- 
verse. Turbulence is characterised by the co-existence of 
motions with many length and time scales described by 
many degrees of freedom. The key to our understand- 
ing of turbulence is to elucidate physics of interactions 
between large scales (eg. large eddies) and small scales 
(eg. turbulent fluctuations), and to develop mathemat- 
ical models that account for the effects of small scales 
without actually solving for them. According to the Lan- 
dau description, superfiuid 4 He consists of the ground 
state and the excitations — quasiparticles drifting on top 
of the ground state. In the language of relativistic quan- 
tum fields this corresponds to vacuum and matter (eg. 
gravity waves interacting with a vacuum). Indeed, there 
is a close relationship between superfiuid hydrodynamics 
and quantum gravity, so that at some level of hierarchy 
of parameters the interactions of the quantum vacuum 
and matter can be described by the defocusing nonlinear 
Schrodinger (NLS) equation [l|. The dynamics of Bosc 
condensates depends on the energy exchange between the 
condensed and non-condensed parts of the gas. Again, 
the NLS equation (reformulated as the Gross-Pitaevskii 
(GP) equation 0) describes equilibrium and dynamical 
properties of BEC as well as the formation of BEC from 
a strongly degenerate gas of weakly interacting bosons 
0, 0] • The formation of the large-scale coherent localized 
ground state (condensate) from a non-equilibrium initial 
state has been studied in a number of papers addressing 
different stages of the formation. Weak turbulence the- 
ory has been used to predict the self-similar evolution of 
the field in the regime of random phases of Fourier ampli- 
tudes [B|, Q , the transition from the regime of weak tur- 
bulence to superfiuid turbulence via states of strong tur- 



bulence in the long- wavelength region of energy space 0] , 
and the final stage, resulting in the formation of a genuine 
condensate [H, 0]. The related question about the effect 
of finite temperature on the BEC dynamics has also been 
addressed recently [Tcj. For instance, it was shown that 
the presence of the thermal cloud in a trapped conden- 
sate creates an effective dissipation that forces a single 
vortex to move away from the centre and disappear [11 1 . 



The problem of the vortex tangle interacting with the 
normal fluid (thermal cloud) is the key question in su- 
perfiuid turbulence. The Landau two-fluid theory of su- 
perfluidity pre-dated the discovery of quantised vortex 
lines and therefore omitted significant dynamical effects. 
This was remedied — in the limit in which the mean spac- 
ing between the vortex lines is small compared with any 
other length scale of interest — by HVBK theory 12, Hf . 
In this limit, the superfiuid vorticity is treated as a con- 
tinuum, but the discrete nature of the vorticity gives rise 
to an extra force on the superfiuid component, arising 
from the tension in the vortex lines. This term is absent 
from the classical Euler equation of motion for an inviscid 
fluid. The vortex lines also create a force of mutual fric- 
tion between superfiuid and normal fluid in addition to 
the mutual friction included by Landau in his equations, 
and represents the effects of collisions of the quasiparti- 
cles with the vortex cores. Such forces were introduced 
into the Landau model in an ad hoc way. This Letter 
is the first attempt to study the effect of these collisions 
quantitatively: we shall find the vortex line decay law 
at non-zero temperature in the context of the defocus- 
ing NLS equation. The NLS equation is a good starting 
point, as the non-dissipative Landau two-fluid model can 
be obtained from the equations of conservation of mass 
and momentum for a one-component barotropic fluid us- 
ing a general expression for the internal energy functional 
of the density [lj] . Through the Madelung transforma- 
tion the NLS equation can be written in that form. Anal- 
ogously, the transport coefficients in the Landau model 
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have been obtained directly from the NLS equation by 
following the Chapman-Engskog expansion [15j. Note 
that the separation of scales needed to carry out the 
derivation of the Landau two-fluid model from the NLS 
equation does not allow the inclusion of vortices as part 
of the ground state. It is natural, therefore, to attempt 
to derive the corresponding effects of the interactions of 
vortices with the quasiparticles directly from the NLS 
equation. 

We consider the normalised defocusing NLS equation 
for the complex function ip 0): 
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The dynamics conserves the total number of parti- 
cles N = J \ip\ 2 dx., and the total energy E = 
J (|Vt/>| 2 + ^l^l 4 ) d*L. We consider the uniform discrete 
system of volume V = TV 3 , which is a periodic box on a 
computational grid with 128 3 discrete points. 

Our goal is to determine the universal decay law for 
the vortex line density in the entire range of tempera- 
tures from to the critical temperature of condensation, 
T\. Our approach consists of three essential steps. We 
aimed to: (1) achieve the thermal equilibrium state for 
the given number of particles and given energy, starting 
from a non-equilibrium stochastic initial condition for the 
wavefunction tp; (2) introduce a vortex ring into this state 
and follow its decay via interactions with non-condensed 
quasiparticles; (3) relate the decay rate to the tempera- 
ture at equilibrium, where we derive the expression for 
the relative temperature, T/T\, as a function of the total 
number density, p = N/V, and the number density of 
the condensate, p . 

We performed large scale numerical simulations of 
Eq. {1} starting from a strongly non-equilibrium initial 
condition Q, where the phases of the complex Fourier am- 
plitudes dk(£) = / ip(x, t)e~ tp x g?x are distributed ran- 
domly at t = 0. Here the momentum p takes quantised 
values p = (27r/7V)n with n = (0, 0, 0), (±1, 0, 0), • • ■. 

This initial state describes a weakly interacting Bose 
gas that is so rapidly cooled below the critical BEC tem- 
perature that the particles remain in a strongly non- 
equilibrium state. The kinetics of the initial weak tur- 
bulent state has been analysed in [f| 0, [Jj| discovering a 
quasiparticle cascade from high energies to low energies 
in the wave number space. The ordering of the system 
and the violation of the assumptions of weak turbulence 
occurs very rapidly in a low-energy part of the spectrum, 
with the formation of a quasi-condensate consisting of a 
tangle of quantised vortices. The vortex tangle decays as 
the system reaches a state of thermal equilibrium with 
some portion (po = |ao| 2 /V) of particles occupying the 
zero momentum state (genuine condensate) and the rest 
of the non-condensed particles being distributed accord- 
ing to the Raylcigh-Jeans equilibrium distribution 17| . 
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where T is the temperature and lob(p) is the Bogoli- 
ubov dispersion relation (see below). An ultraviolet cut- 
off for this distribution appears naturally through the 
spatial discretization of the NLS equation. The numer- 
ical scheme consists of fourth-order finite difference dis- 
cretization in space and fourth-order Runge-Kutta in 
time, so it is globally fourth-order accurate. This scheme 
corresponds to the Hamiltonian system in the discrete 
variables ipjkn, such that 
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j, k,n = 1, ...,7V , (3) 
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where 



H = Efe-fe^ " 1*1 + f^n] + i|</W| 4 )(4) 
jkn 

with * 2 = 1pj+2,k,n + tpj-2,k,n + tpj,k+2,n + 1pj,k-2,n + 
i'j,k,n+2 + ipj,k,n-2 and *i = ipj+i^.n + ipj-i,k, n + 

1pj,k+l,n + Ipj ,k— l,n + 1pj,k,n+l + 1pj,k,n-l- 

The thermodynamic description of the condensation 
process has been obtained in [9( by adapting the Bogoli- 
ubov theory of a weakly interacting Bose gas [l8| to the 
classical system ([!]). We follow the same basic idea to de- 
rive expressions for the energy and non-condensed part of 
the discretised energy Q written in terms of the Fourier 
amplitudes a p as 
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where 5 P is the Kronecker delta symbol and 



*a(P) = oE sin2 (P'/ 2 )( 7 " cosfe))- 



(5) 
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The Bogoliubov transformation b p = u p a p — v p a*, 



such that Up = 1/y 1 — Q P and v p = Q p / yjl — Q p with 
Q p = [—K2 — 2po+ub(p)]/po diagonalises the term in (JSJ), 
which is quadratic in ao, to ^2 p 'uJb(p) b p b p , where Y^' p 

excludes the p = mode. Here wb{p) = \J K% + ^Po^2 
is the Bogoliubov-type dispersion relation. 

Using the equilibrium distribution of the uncondensed 
particles j2]) the non-condensed number density can then 
be expressed in terms of the basis used in this diagonal- 
isation as 



modified by the presence of nonlinear interactions with 



T^K 2 (p)+p 



(7) 
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The discretised energy density H/V in the new basis 
takes the form 



H 



i[p 2 + (p-po) 2 ]4£'i. 
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The Eqs. ©-© are analogous to Eqs. (8)~(9) of [§] but 
modified for the discrete Hamiltonian discretization 
Given the energy density. H/V, and the total number 
density, p, one can determine the temperature, T, at 
equilibrium and the number density of the condensed 
particles, po, from Eqs. j7|) and (J8]). The condensate 
fraction po/p as a function of the energy density H/V 
is shown in FIG.l. This figure can be compared with 
FIG. 2 of [j| for the spectral representation of the total 
energy. The analytical formulae ©-I© predict the sub- 
critical behaviour of condensation, whereas the numerics 
does not support this conclusion, as shown in the insert 
of FIG.l. We use a linear approximation for small po to 
determine the critical maximum energy for condensation 
as shown in the insert. This energy is then used to deter- 
mine the critical temperature for condensation T\ (= T 
for min _ff/T^ for which po = 0) from (J?])-©. We found 
a phcnomcnological formula that determines T/T\ as a 
function of po and p as 

l-'-O-^f-^?)'. < 9 > 

where a is the only fitting parameter that we found as 
a = 0.227538. The insert in FIG.l shows the graph of 
T/T\ as a function of po/p for p = 1/2. Eq. ([9]) gives an 
excellent fit to the values computed from (J?])-© across 
all the values of po and p. 

In order to analyze the decay of the vortex line length 
at non-zero temperatures, we insert a vortex ring into a 
state of thermal equilibrium and follow its decay due to 
the interactions with the non-condensed particles. The 
condensate healing length, which determines the size 
of the vortex core, is calculated based on the density 
of the condensate, and in our non-dimensional units is 
£ = I n neann g lengths, the radius of the ring 

is set to Rq = 10. The new initial state is ip v (t = 
0) = ip cq * ^Vortex, where ^eq is the equilibrium state 
and ^Vortex is a wavefunction of the vortex ring [l^| . The 
vortex line length, L, is calculated as a function of time 
with high frequencies being filtered out from the field ijj, 
according to a p = a p * max( y/ 1 — p 2 /p 2 , 0) , where the 
cut-off wavenumber is chosen as p c = 1Q(2-k/N) [2Cj| . 
The first important conclusion of our numerical simula- 
tions is that at all temperatures, the square of the vortex 
line length decays linearly with time, 



(10) 



where 7 does not depend on t. FIG. [2] shows this depen- 
dence for various temperatures. The actual isosurfaccs 
of the decaying vortex line are shown in the inserts. 



FIG. 1: (colour online) Condensate fraction, po/p, as a func- 
tion of the energy density as obtained from the numerical 
simulations (points) and from the analytical expressions (J7J- 
© (solid line). The inserts show (a) the plot of T/T x as 
a function of po/p, obtained using Eqs. (|7]l-(|8)l or Eq. @ 
and (b) subcritical condensation predicted by Eqs. 0-© 
(black line), the linear approximation used to obtain the crit- 
ical temperature of condensation (gray (red) line), and nu- 
merical calculations (blue dots). The total number density is 
p=l/2. 
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FIG. 2: (colour online) The decay of the square of the vortex 
line length as a function of time for various T indicated next 
to the graphs. The fit to the linear function is shown by 
the gray (red) lines. The inserts show isosurface plots of the 
vortex line (for filtered fields ip; see text) for T = 0.52T\ at 
time=130 (left) and time=1300 (right); between these two 
times the vortex line length is reduced by a factor of 2. The 
perturbations to the vortex line due to collisions with non- 
condensed particles are clearly seen on the left insert. These 
collisions generate Kelvin waves that also radiate energy to 
sound. The total number density is p — 1/2. 




100 200 300 400 500 600 4-\ 



time 



This result agrees with predictions of the HVBK theory 
for supcrfluid helium 12[ according to which the funda- 
mental equation of the motion of a vortex line, vl, is 
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given by (see also page 90, Eq. (3.17) of [211 ] ) 

v L = v s i + as'x (v n - v s i)-aVx[s'x(v n - v s i)], (11) 

where v s i is the local superfluid velocity that consists of 
the ambient superfluid flow velocity and the self-induced 
vortex velocity Ui, v n is the normal fluid velocity, s is 
a position vector of a point on the vortex and s' is the 
unit tangent at that point. Mutual friction parameters a 
and a' are ad hoc coefficients in the HVBK theory that 
are functions of p n , /?, and T only. Eq. is a gen- 
eral and universal equation used to follow the evolution 
of three-dimensional vortex motion in an arbitrary flow. 
When formulated for a single vortex ring Eq. [TT] reads 
dR/dt = -am, where m = n[log(8R/^) - S + l]/(47i\R) 
and S is the vortex core parameter. For the GP vortices 
5 w 0.38 Q. In dimensionless units used in our paper 
Ui = [log(8i?) — S + l]/R. After integration of the equa- 
tion for R we get at = (R 2 - R 2 )/[2(\og(8R) + 5 - 1)], 
where R is the mean radius of the ring. When this is com- 
pared with (TIT)]) we get the following relationship between 
7 and a: 7 = 87r 2 (log(8i?) + 6 — l)a. From our numer- 
ics we obtained a general result valid across all ranges of 
temperatures and total densities: 7 w Kp(T/T\) 2 , where 
K fn 68. Note that for a GP condensate T/T x ~ p n /p 
to the first order (see insert (a) of FIG.l), so alterna- 
tively, we can write 7 k, K\p n (T /T\) FIG. [3] shows the 
comparison of the numerically calculated 7/p and the 
quadratic fit K(T/T\) 2 . Thus, we found that the mutual 
friction coefficient in condensate superfluids is given by 
anK 2Pn (T/T x ). 

The existence of the transverse force on superfluid vor- 
tices which is parametrised by the parameter a 1 has been 
a subject of much debate in mid-1990s, when calcula- 
tions of the classical Magnus force applied to superfluid 



vortices have been offered and argued about [22|. The 
criticism is based on the observation that the classical 
hydrodynamic equations are inapplicable in the vortex 
core. Whether or not the details of the non-classical vor- 
tex dynamics are crucial to the existence of the trans- 
verse force is still an open question. The estimate of a' 
can be obtained from our numerical procedure as follow- 
ing. Eq. (fTTj) written for a distance travelled by a single 
vortex ring takes form (see Eq. (3.53) on page 107 of [2l| ) 
dz/dt = (1 — a')itj. We compared the distances travelled 
by a vortex ring at various temperatures obtained nu- 
merically with the distances travelled by a vortex ring in 
the absence of the transverse force according to the ana- 
lytical formula dz/dt = Ui, where Ui — Ui(R(t)) and R(t) 
varies with time according to (TIT))) . The insert of FIG(3] 
shows these distances for T/T\ = 0.27. Our calculations 
fail to detect any significant presence of the transverse 
force for any temperature considered: the deviation from 
the analytical curve is insignificant within the accuracy 
of (HUl). We plan to perform a more thorough analyti- 
cal and numerical study of transverse force from a single 



FIG. 3: (Color online) Values of 7/p as a function of tempera- 
ture T/Tx for various values of the total number density p de- 
picted in various shades of gray (in various colours) : p = 1/2 
(dark (red)), p = 1/4 (light (green)) and p — 3/4 (medium 
(blue)). The plot of the quadratic fit j/p = 68(T/T A ) 2 is 
given by the dashed line. The relative temperature is calcu- 
lated using Eq. ([5}. The result is not sensitive to whether 
we use the values of p and po that correspond to the state of 
thermodynamical equilibrium before the introduction of the 
vortex ring or after the vortex ring disappears and the sys- 
tem equilibrates. The insert shows the distance travelled by 
a vortex ring as a function of time for T/Tx = 0.27 (red dots 
- distances calculated using dz/dt — m, black line using nu- 
merics). Curves depart when the vortex ring becomes small 
in radius and the analytical formula is no longer accurate ap- 
proximation of the vortex velocity. 
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phonon acting on a single vortex in context of the GP 
model in future. 

In summary, we considered the effect of temperature 
on the decay of vortex line density via interactions with 
non-condensed particles in the context of the defocus- 
ing NLS equation. We obtained a simple expression for 
the temperature at equilibrium as a function of the to- 
tal number density and the number density of the con- 
densed particles. Depending on these two parameters, 
a vortex ring introduced into the condensate shows dif- 
ferent decay rates with time. We identified this decay 
law as linear for the square of the vortex line length and 
showed the universal dependence of the decay rate on 
temperature and total density. It has been suggested 
that the emission of sound by vortex reconncctions and 
vortex motion is the only active dissipation mechanism 
responsible for the decay of superfluid turbulence. The 
decay of superfluid turbulence via Kelvin wave radiation 
and vortex reconnections was studied in the framework 
of the GP equation [23| at near zero temperature, via 
collision of two vortex rings, and confirmed that in the 
Kelvin wave cascade, where energy is transferred to much 
shorter wavelengths with a cut-off below a critical wave- 
length, the vortex line density can be described by the 
famous Vinen equatio n | 24|| d(L/V)/dt = — x(L/V) 2 . It 
has also been shown [25j that the presence of localized 
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finite amplitude sound waves greatly enhances the dissi- 
pation of the vortex tangle, essentially changing the de- 
cay law to exponential decay. This Letter complements 
the existing Kelvin wave cascade scenario by consider- 
ing an opposite limit when there are no reconnections, 
and the decay mechanism depends only on the energy 
exchange with non-condensed particles. This mechanism 
exceeds the energy transfer via the Kelvin wave cascade. 
Finally, we related our results about a single vortex ring 
to the mutual friction coefficients in the general equation 
of vortex motion. 
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